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ABSTRACT 

We report on spectral and timing properties of the magnetar CXOU J164710. 2— 455216 in the 
massive star cluster Wester lund 1. Using 11 archival observations obtained with Chandra and XMM- 
Newton over approximately 1000 days after the source's 2006 outburst, we study the flux and spectral 
evolution of the source. We show that the hardness of the source, as quantified by hardness ra- 
tio, blackbody temperature or power-law photon index, shows a clear correlation with the 2-10 keV 
absorption-corrected flux and that the power-law component fiux decayed faster than the blackbody 
component for the first ^^100 days. We also measure the timing properties of the source by an- 
alyzing data spanning approximately 2500 days. The measured period and period derivative are 
10.610644(17) s (MJD 53999.06) and < 4 x lO^^^ g g-i (99% confidence) which imply that the spin- 
inferred dipolar magnetic field of the source is less than 7 x 10^^ G. This is significantly smaller than 
was suggested previously. We find evidence for a second fiux increase, suggesting a second outburst 
between MJDs 55068 and 55832. Finally, based on a crustal cooling model, we find that the source's 
cooling curve can be reproduced if we assume that the energy was deposited in the outer crust and 
that the temperature profile of the star right after the 2006 outburst was relatively independent of 
density. 

Subject headings: pulsars: individual (CXOU J164710. 2— 455216) - stars: magnetars - stars: neutron 
- X-rays: bursts 



1. INTRODUCTION 

Neutron stars with ultrastrong magnetic fields, so- 
called magnetars, can emit radiation which is or- 
ders of magnitude stronger than their rotation power. 
These objects are prone to X-ray outbursts in which 
the flux r ises rapidly then decays over days to 

• months (seelWoods Thompson 2006; Mereghctti 2003; 

I iRea fc Espositoll201lL for reviews). Such outbursts are 
often accompanied by short X-ray and soft gamma-ray 
bursts. The origins of outbursts are still debated. It has 

I been suggested that a sudden or gradual internal heat 
release possibly caused by crustal cracking, heats the 
star from within, and induces observed temperature and 

] luminosity increases ()Thompson fc Duncan! [19951 [1991 
iPerna &: Pons! [20T1I) . It has also been suggested that 
crustal moti on can generate twists in the external mag- 
netic fields ([Thompson et al.l l2002t ) . The twisted mag- 
netic fields induce currents in the magnetosphere, which 
return to the stellar surface and heat it. 

The flux relaxation following an outburst can 
be explained as passive cooling of the hot star 
([Lvubarskv et al.|[2Q02[: [Pons fc Reall20ia see also Gum- 
ming et al. 2012, in preparation), or as "u ntwisting" 
of the external magnetic fields (Be loborodovl 2009). In 
either case, fiu x and spectral hardnes s are expected to 
be correlated ([Thompson et al.l l2002t iLvutikov I [20031: 
lOzel fc Giiveii I2007D . The correlation has been ob- 
served in many magnetar s (e.g., S GR 180 6—20, IE 
2259+586, IE 1547-5408: Woods et al. 2007; IZhu et al.l 
I2008t iScholz fc Kaspill201l[ ) afthough it is not as clearly 

^ Lome Trottier Chair; Canada Research Chair 



seen in some sources (e.g., SGR 19 00-H4, SGR 1627-41; 
ITieneo et all [20071 : 1 An et al.l[20Tl . 

The high magnetic fields of magnetars have been in- 
ferred both from spin-down rates, assuming the stan- 
dard dipole braking relation B ee 3.2 x 10^^ {PP)^/'^ G 
([Manchester fc Tavloil[T977l ) . as well as from indirect ar- 
guments regarding rad iative behavior and field decay 
(iThompson fc DuncanI [19951 [19960 . Indeed the typical 
magnetic field for an object classified as a magnetar on 
the basis of radiative behavior is 10"'^* — 10"'^^ G[3 However 
recently a small handful of objects have been reported 
to have fields below this range, overlapping with those 
of apparently ordinary radio pulsars (S GR 0418-1-5729, 
Swift J1822 .3-1606; React al. 2010; L ivingstone et al.l 
[Mil: iRea et all l20ia IScholz et al.l 12012). Although 
the i3-distribution of radiatively classified magnetars 
still re mains significantly highe r than that of radio 
pulsars ([Livingstone et al] 1201 ID . the apparently low- 
B magnetars are puzzling and are suggestive of higher 
order multipolar structure or strong internal toroidal 
fields, for which the dipole component is deceptively low 
([Turolla et al.l[20H . 

CXOU J164710.2-455216 was discovered on 1998 June 
15 with the Chandra X-Ray Observatory ([Muno et al.l 
[200I) . It is located at R.A. = 16M7"'10".2 0, Deck = 
-45°52'17".05 rj2000.0. [Skinner et al.l[2006[) and is esti- 
mated to be 2. 5 — 5 kpc away, in the massive star cluster 
Westerlund 1 ([Clark et al.ll2005D . A short 20- ms burst 
was detected with the Swift Burst Alert Telescope (BAT) 



See the onhne magnetar catalog for a 

compilation of known magnetar properties, 

[http : / /www. physics. mcgill.ca/~pulsar/magnetar/main. html 



Table 1 

Summary of observations and spectral fit results 



# 


Date 


Observatory 


Mode 


ID 


kT 


r 


rSrS radius 


Flux 


T5T trim^r^ 

r L c mx 


AT" 




M JD 








(kcV) 




(km) 






(s) 


1 


53513.0 


Chandra 


TE 


6283 


0.51(2) 




0.50(5) 


0.030(2) 




3.24 


2 


53539.9 


Chandra 


TE 


5411 


0.53(14) 


3.33(55) 


0.25(13) 


0.028(11) 


0.02(1) 


3.24 


3 


53995.1 


XMM 


FW/FW 


0404340101 


0.59(6) 


3.86(22) 


0.18(4) 


0.025(4) 


0.016(3) 


2.6/0.073 


4 


54000.7 


XMM 


FW/FW 


0311792001 


0.70(1) 


2.90(6) 


1.44(5) 


3.00(7) 


1.62(6) 


2.6/0.073 


5 


54005.4 


Chandra 


cc 


6724 


0.60(1) 


2.46(13) 


2.22(9) 


2.52(9) 


0.99(7) 


0.00285 


6 


54010.1 


Chandra 


cc 


6725 


0.60(1) 


2.69(12) 


2.08(8) 


2.08(7 ) 


0.74(6) 


0.00285 


7 


54017.4 


Chandra 


cc 


6726 


0.61(1) 


2.92(13) 


1.97(6) 


1.84(6) 


0.53(5) 


0.00285 


8 


040OD.4 


Chandra 




8455 


0.58(1) 


O "70 /" 1 '7^ 

2. /o(l / ) 


1.97(9) 


1.42(6) 


0.41(5) 


0.00250 


9 


54133.9 


Chandra 


cc 


8506 


0.56(1) 


2.86(18) 


1.83(9) 


0.95(4) 


0.29(3) 


0.00285 


10 


54148.5 


XMM 


SW/LW 


0410580601 


0.58(1) 


3.20(11) 


1.39(5) 


0.78(2) 


0.27(2) 


0.3/0.048 


11 


54331.6 


XMM 


SW/LW 


0505290201 


0.58(1) 


3.42(10) 


0.99(4) 


0.42(2) 


0.17(1) 


0.3/0.048 


12 


54511.5 


XMM 


SW/LW 


0505290301 


0.54(2) 


3.40(14) 


0.93(7) 


0.26(1) 


0.11(1) 


0.3/0.048 


13 


54698.7 


XMM 


SW/LW 


0555350101 


0.57(2) 


3.76(15) 


0.62(4) 


0.16(1) 


0.068(8) 


0.3/0.048 


14 


55067.6 


XMM 


SW/LW 


0604380101 


0.53(2) 


3.68(15) 


0.56(4) 


0.087(5) 


0.036(4) 


0.3/0.048 


15 


55832.0 


XMM 


SW/LW 


0679380501 


0.73(1) 


2.89(11) 


0.67(3) 


0.54(2) 


0.18(2) 


0.3/0.048 


16 


55857.8 


Chandra 


TE subarray 


14360 


0.62(2) 


1.50(47) 


0.87(8) 


0.44(4) 


0.18(2) 


0.44 



Notes. Tfie outburst occurred on 53999.05659 (MJD). Obs. 3 (MJD 53995.1) was used to set the quiescent level for tlie flux evolution 
(Fig.[TJ. TVh viSts obtained by simultaneously fitting Obs. 4-14, and was fixed for fitting the others. Fits are conducted in the 0.5-10 keV 
band, and uncertainties are at the 1-a confidence level. 

^ MOSl,2/PN for the JfMM-Afe-iuion observations. TE: Timed Exposure, CC: Continuous clocking, FW: Full Window, LW: Large Window, 
SW: Small Window. 

^ Blackbody radius. Results of tbabs (bbodyrad+pow) fitting in XSPEC for an assumed distance of 5 kpc. 



Absorption-corrected fiux in the 2-10 keV band in units of 10 erg cm 
Absorption-corrected power-law flux in the 2-10 kcV band in units of 10^ 
° Time resolution. 



erg cm 



on 2006 Sept. 21 (MJD 53999.05659. lKrimm et aL|[200l . 

The source was subsequently observed with multiple X- 
ray observatories i ncluding Sw ift, XMM-N ewton, Chan- 
dra an d Suzaku. Ilsrael et "all (^007) and iWoods et al.l 
(j2011h reported on the data obtained after the 2006 out- 
burst out to approximately 100-200 days and measured 
spectral and timing properties, where they concluded 
that there is no clear evidence of significant s pectral vari- 
abilit y. The pulsar's spin period is ~ 10.6 s (jMuno et alJ 
[200i l and the spin-down rate has been reported to be 
0.8 - 1.3 X 10-12 s s-i (|Israel et al.ll2007t IWoods"eraII 
1201 Ih . These previously reported spin parameters im- 
ply an inferred surface dipolar magnetic-field strength of 
B = 0.95- 1.2 X 10^^ G. 

Here, we report on the analysis of Chandra and XMM- 
Newton archival data which cover a longer time period 
(^2500 days). We study the long-term spectral evolution 
of the source after its outburst and report on the pulsar's 
timing behavior during this period. 

In Section O we summarize the observations we used 
for our study. We then describe our data analysis and 
results in Section [3l discuss the results in Section IH and 
finally present our conclusions in Section [5l 

2. OBSERVATIONS 

Table [1] presents the data we used for our analy- 
sis. In total, we analyzed 16 observations obtained 
with the Chandra and XMM-Newton telescopes over 
the course of '^2500 days. For the Chandra observa- 
tions, we reprocessed the standard pipeline output using 
chandrajrepro of CIAO 4.43 along with CALDB 4.4.7 to 
produce the "level 2" event lists using the newest soft- 
ware and calibration updates^ For the XMM-Newton 
data, we processed the Observation Data Files (ODF) 
with epproc and emproc and then applied the standard 

http:/ /cxc. harvard. edu/ciao4.4/index. html 
* http:/ /cxc. harvard.edu/ciao/threads/createL2/ 



filtering procedure (e.g., flare rejection and pattern selec- 
tion) of Science Analysis System (SAS) version 11.0.00 
Some observat ions in the Table were already analyzed by 
other authors (iMuno et al.l l2006l 120071: Ilsrael et al.ll2007i: 
IWoods et al.ll201lD . However, we re-analyzed those using 
the above procedure for consistency. 

3. DATA ANALYSIS AND RESULTS 
3.1. Imaging Analysis 

We detected the source with wavdetect for the Chan- 
dra TE mode observations and with edetect_chain for 
the XMM-Newton observations. The source positions we 
found are all consistent with one another and with the 
known source location within 90% uncertainties of the 
position measurements (0".6 for Chandr^ and 2" for 
XMM-Newtoi^. 

The relatively large point-spread function of XMM- 
Newton was a concern since the source is located in a 
star cluster, and any other object within ~ 30" could 
in principle contaminate the source spectrum. With 
two Chandra TE observations (IDs 6283 and 14360), we 
searched for X-ray sources within a radius of 30" cen- 
tered at CXOU J164710.2-455216 but found none. We 
also checked whether our source is consistent with being 
a point source using the eradial tool for XMM-Newton 
observations and Chandra Ray Tracei0 and the MARjfl 
tools for the Chandra observations. The data were con- 
sistent with the target being a point source in all the 
observations except for the first XMM-Newton observa- 
tion (ID: 0311792001), where the radial profile of the 
source was distorted due to pile-up. Therefore, we con- 
clude that there is no significant contamination from un- 

http://xmm.esac.esa.int/sas/ 
^ http://cxc.harvard.edu/cal/ASPECT/celmon/ 
''' http:/ /xmm2.esac.esa.int/external/xmm_sw_cal/calib 
* http://cxc. harvard . edu/chart / 
9 http://space.mit.edu/CXC/MARX 
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resolved sources in the XMM-Newton observations and 
that the source exhibits no detectable extended emission. 

3.2. Spectral Analysis 

We used all the data listed in Table [1] for our spectral 
analysis. For the Chandra observations, we extracted the 
source events using a box of dimension 2" along the 1-D 
events and 5" in the orthogonal direction for the CC- 
mode observations, and a circle with radius 2" for the 
TE-mode observations. Backgrounds were obtained in 
two rectangular regions with size of 10" x 5" from each 
side of the source region and an annular region with radii 
of 5" and 10" centered at the source for the CC-mode 
and the TE-mode observations, respectively. Then we 
produced spectra using the specextract tool of CIAO 
4.4 with CALDB 4.4.7 and grouped them to have a min- 
imum of 20 counts per bin for further analysis. 

For the XMM-Newton data, we extracted source spec- 
tra from circular regions having radius of 16" and back- 
ground spectra from source-free regions on the same chip. 
Corresponding response files were produced using the 
rmf gen and the arf gen tasks of SAS 11.0.0. Each spec- 
trum was then grouped to have a minimum of 20 counts 
per bin. 

The first XMM-Newton observation was mildly piled 
up, so we excluded the central ~5" in each of the PN and 
the two MOS detectors to minimize the pile-up effect. 
After removing the central regions, we checked if pile-up 
was still significant using the epatplot tool of SAS. With 
the removal, the measured event pattern distributions 
showed good agreement with the expected ones. We also 
checked if the spectral parameters changed significantly 
when we removed central regions of different sizes. The 
power-law index changed smoothly by approximately 3% 
when we varied the removal radius from 0" to 9" for the 
PN data, but it changed abruptly (softened by 10%) for 
removal regions between 0" and 3". 5 and then stayed 
constant for the MOS data. Therefore, we ignored the 
central 5" for both the PN and the MOS data to ensure 
that pile-up would not distort the spectrum of the source 
in this one observation. 

We fit the 0.5-10 keV spectra with three models: an 
absorbed power law plus blackbody, an absorbed double 
blackbody and an absorbed blackbody (tbabs* (power 
+ bbody), tbabs* (bbody + bbody) and tbabs*bbody 
in XSPEC 12.7.0)E3 We fit aU the data in which the 
source flux decreased monotonically after the 2006 out- 
burst (Obs. 4-14; see Table [T]) simultaneously with a 
common hydrogen column density (iVn)- 

The single blackbody fit {x^/DoF ^ 9117.42/7062) 
and double blackbody fit {x^/DoF = 7353.82/7040) 
were not acceptable as they showed systematic trends 
in the low energy (< 1 keV) and the high energy 
bands (> 8 keV). However, the blackbody plus power- 
law fit described the observed spectra well [y^ / DoF = 
6988.04/7040). The iVn we obtained from the fit is 
2.39(5) X 10^^ cm~^, which is different from the val- 
ues that [Muno et al.l H2007I . 1.44 xlO^^ cm'^) and 
[Israel et all (|2007l 1.9 xlO^^ cm'^) use d. However 



it is consistent with that obtained bv iWoods et al.l 
(|201H) and that reported by iSkinner et all (|2006l 2.~ 
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Figure 1. Flux evolution of CXOU J164710.2-455216 after the 
2006 outburst and crustal cooling model fits. Symbols are 2-10 keV 
absorption-corrected flux evolutions for total (diamonds), power- 
law (triangles) and blackbody (squares) components, and lines are 
quiescent flux (dot-dashed, Obs. 3 in Table[l} an d cru stal cooling 
models for decay of the source flux (see Section 14.21 for details). 



Solid line: Tc = 1.1 X 10^ K, B = 8 X lO" G, Qi, 



5, initial 



http: / /heasarc. gsfc.nasa.gov/xanadu / xspec / 



temperature changes from 1.37 X 10 K at 9 X 10 g cm to 
4.5 X 10* K at 5 X lO" g cm'^. Dotted line: Tc = 9 X lO'' K, B = 
2 X 10^* G, Qimp = 1, initial temperature changes from 1.3 X 10^ K 
at p = 1 X 10^ g cm^^ to 4 X 10* K at p = 5 X lO^^ g cm"^. 
Dashed line: Tc = 5.7 X 10^ K, B = 2 X lO^^ G, Qimp = 5, initial 
temperature is T = 1.5 X 10^ K for 1.5 X 10^ g cm"^ < p < 
1 X 10^^ g cm~^. Other parameters were held fixed: M = 1.4Mq, 
_R = 12 km. Note that the uncertainty for each data point is smaller 
than the symbol size. 



3.0x10^ ^ cm~^) on t he ba sis of cluster extinction esti- 
mates oflClark et al.l (|2OO50 and the empirical relation of 
iGorensteinI (|1975| ). We then fixed iVn and tried to fit the 
spectra of the other, fainter observations with the same 
models. The results are summarized in Table [TJ 

We assumed that the source was in (or near) quies- 
cence for the first three observations (Obs. 1-3). The 
source spectrum was well fitted with a single blackbody 
for Obs. 1 without requiring an additional component, 
although a power law was not ruled out unambiguously. 
We fit Obs. 2 and 3 with a power law; a single blackbody 
fit was unacceptable. However, adding another compo- 
nent (blackbody) significantly improved the fit (F-test 
probabihties of 0.048 and 0.00063 for Obs. 2 and 3, re- 
spectively). 

Another flux increase was observed in Obs. 15, which 
might be due to a second outburst between MJDs 55068 
and 55832 (see Fig. [T]). We were not able to fit the 
spectrum for Obs. 15 or 16 with single-component mod- 
els. For Obs. 15, a blackbody plus power law pro- 
vided a better fit (p = 0.080) than a double blackbody 
{p = 0.033 for kTi = 0.36(2) keV, /cTa = 0.82(2) keV). 
For Obs. 16, a double blackbody model fit the data bet- 
ter with kTi = 0.58(3) keV and fcTz = 1.33(35) keV 
although a blackbody plus power law also provided an 
acceptable fit. Since it seems unlikely that the black- 
body temperature increased from 0.80 keV to 1.33 keV 
without any accompanying flux increase, we report the 
blackbody plus power-law spectrum for Obs. 15 and 16. 
The results of all our spectral fits are reported in Table [TJ 

We also fitted the spectra for Obs. 3-14 with 
an absorbed resonant cyclotron scattering model, 
tbabs* (atable{RCS. mod}) in XSPEC (,Rea et al..,200&l . 
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Figure 2. The time evolution of various spectral properties. The top axis shows the 2-10 keV absorption-corrected flux interpolated with 
the double-exponential relaxation fit to the flux evolution (see section 13.31 1 . The conversion between time and flux was obtained from the 
double exponential fit result. Clear trends in time (flux) are visible, a) Hardness ratio (i<4_io kov/^2-4 kcv)- f') Power-law photon index, 
c) Blackbody temperature, d) Blackbody radius assuming a distance of 5 kpc. 



From the fit, we obtained a smaller Nn (2.06 ± 0.02 x 
10^^ cm~^), as expected because of the strong cutoff in 
the power-law spectrum at low energies. We were able to 
measure 2-10 keV absorption-corrected fluxes using the 
model, and they agreed well with those of the blackbody 
plus power-law fits described above. However, the RCS 
model parameters were not well constrained. Therefore, 
we mainly use the results of the blackbody plus power- 
law fit for discussion below. 



3.3. Spectral Evolution 

We plot the time evolution of the 2-10 keV absorption- 
corrected fluxes after the 2006 outburst in Figure [TJ 
We have attempted to fit the data out to ^1000 days 
with a double exponential, F{t) — F^e"'*"*")/"^! -I- 
i^2e~*-*"*°-'^^^ + Fq or a broken power law, F(t) ~ 
Flit - t^Y^ + Fq for to < t < Tbreak and F{t) = 
F2(t - to)"' + Fq for Tbreak < t, whcrc Fi and Fa are 
determined so that F{t) is continuous ai t = Tbreak- In 
fitting, we set t^ to be the onset of the outburst (MJD 
53999.05659) and fixed the quiescent flux [Fq) to be the 
measured value of 2.5 x 10~^^ erg cm~^ s~^ (Obs. 3 in 



Table [J). 

Although the fit results with the simple analytical 
models were statistically unacceptable, the double ex- 
ponential function roughly described the relaxation with 
decay times of 18 ±2 and 325 ±10 days (xVdof=33.5/7), 
as did the broken power-law fit with decay indices of 
-0.22 ± 0.01 and -1.16 ± 0.02, and Tbreak of 106 ± 5 
days (x^/dof=64.8/7). Note that we include only statis- 
tical uncertainties in fitting, although some systematic 
for example, the sudden decrease in flux at ^150 days in 
Figure [1] may indicate a cross-calibration issue between 
Chandra and XMM-Newton. Indeed, adding a 6% sys- 
tematic uncertainty to the fluxes made the double expo- 
nential fit acceptable (decay times of 19 ± 4 and 331 ± 15 
days, xVdof=12.3/7). 

We also plot the time evolutions of the spectral proper- 
ties for the same data in Figure [2] in order to look for any 
correlation between spect ral har dness and flux, as seen 
in other magnetars (e.g.. |Rea et al. 200 S l\yoods et al 
2007t iCampan a et~all 120071: ifam et al. .20081 : iZhu et a! 



20081 iScholz fc Kaspill2011[ ). Clear trends are visible in 
the hardness ratio, power-law photon index and black- 
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Table 2 

Timing properties of CXOU J164710.2-455216 





Period 


Period Derivative 


Epocli 
(MJD) 


B° Comment 
(10" G) 


Israel ct al. (2007) 


10.6106549(2) 


9.2(4) X lO"^-* 


53999.0 


1.00 


Woods et al. (■2011') 


10.6106567(1) 
10.6106567(2) 


8.3(2) X 10"^^ 
1.3(1) X 10"^^ 


54008.0 


0.95 

1.19 witli cubie term 


This work 


10.610644(17) 


< 4 X 10""'' 


53999.1 


<0.7'' 



Notes. ^ Spin-inferred dipolar magnetic field. 
90% upper limit. 

body temperature; the spectrum clearly softened as the 
flux decreased. 

For the first ^100 days, we were able to re produce 
the same trends as those of I Woods et al.l (|2011[ ) for the 
spectral parameters, though with small offsets within un- 
certainties in the power-law index and the blackbody 
temperature. The offsets could be due to the fact that 
they used a different CALDB or analysis software version 
and/or included a Suzaku observation. 

3.4. Timing Analysis 

Timing studies have b een done by [Israel et al.l (|2007[ ) 
and Woods et al.l (j2011l ) using phase-coherent analyses 
with data obtained for approximately 100-200 days after 
the 2006 outburst. Since the separation between obser- 
vations 10 and 11 was large, phase connection was lost. 
Instead, we searched for the best pulse period in each ob- 
servation to find the ave rage spin-down rate. To do this, 
we employed the H-test (jde Jager et al.llT989l ) because it 
requires no advance knowledge of the pulse profile at any 
epoch, an d because this source's pulse profile i s known 
to change ijlsrael et al.l 120071: iWoods et al.ll201lD . 

The pulse period of the source is ~10.6 s. The time 
resolutions of the observations are all smaller than ~5 s 
so could be used in our analysis (see Table [T|). Thus, the 
baseline useful for timing spans approximately 2500 days. 
We extracted source events from the regions described in 
Section [321 Then we selected events having energies in 
the range 0.5-8 keV for the analysis. Each observation 
was then corrected to the barycenter using the axbary 
tool of CIAO for the Chandra data and the barycen tool 
of SAS for the XMM-Newton data. 

For each observation, we calculated t he sum of the 
Fourier powers of n harmonics, Z^{P) (jBuccheri et al.l 
|l983) for n=l-20 at test periods in the range 10.6050- 
10.6160 s with step size AP = 10~^ s a nd found n and 
P wh ich maximized H = Z^^ (P) — An + 4 ()de Jager et al.l 
119891) for each observation. We folded the events at the 
best period to obtain a pulse profile. We used simula- 
tions to determine the uncertainty on the best period. 
To do this, we used the measured pulse profiles (which 
had 40 bins) to produce fake event lists containing the 
periodicities. Fake event lists were made by first taking 
an observed pulse profile, and adding Poisson noise as 
well as a random phase offset to it. Then, for each bin 
in the profile, we took the total number of counts, as- 
signed each one a random arrival time in that phase bin 
and rebinned the arrival times to the time resolution of 
the detector. In this way, we constructed a simulated 
event list. For each such list, we then calculated the H 
statistics for each P in the range given above, and for all 
n's in the range of 1-20. We then found P and n that 
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Figure 3. The evolution of pulse perio d afte r the 2006 outb urst. 
Periods are measured using the H-test Ide Jag er et al.lll989l ) and 
uncertainties are estimated using simulations (see text). The solid 
line is the best fit to the data, the triple dot-dashed line is the 
90% upper limit for the best fit, and the dotted and the dashed 
lines are the expected evolution of the period in time for the pre- 
viously repo rted period derivative wit h 90% confidence interval by 
I Israel et"aJ] f2007) and IWoods et al.l l[201f ) with no cubic term, 
re spectively. Also plotted in the dot-dashed line is the solution 
of lWoods et al.l (i201fi ) assuming glitch recovery following the 2006 
outburst. 

maximized H, and took the standard deviation of the 
periods as the period uncertainty. The uncertainties we 
obtained in thi s way agree w i th tho se calculated using 
the method of iRansom et al.l (|2002D within a factor of 
two, but typically are slightly larger. To run the simula- 
tions for observations with coarse time resolution where 
the measured pulse profiles were not well determined, we 
assumed the XMM FN profile. 

In Figure [21 we show the best period we obtained for 
each observation. For the XMM-Newton observations, we 
took the variance-weighted average period for measure- 
ments done with the three instruments (MOSl, M0S2 
and PN). We then fit the periods with a linear function 

{P{t) = Po-\-Pt) to obtain the best period at the time of 
outburst, Pq and the spin-down rate. The fit results were 
Po = 10.610644 ± 0.000017 s and P = 0.7 ± 2.4 x 10"" 
s s"^. The period derivative we measured is consistent 
with zero, so we report the 90% upper limi t, 4 x 10~^^ 
s s"^. This is s ignificantly smalle r than what [Israel et al.l 
(|2007[ ) and W oods et all (|201lD reported (see Table [1 
and Section [4?3l for discussion). 

We also searched for short-term aperiodic variability of 
the source in the 0.5-8 keV band. For the XMM-Newton 
data, we produced light curves using evselect and cor- 
rected them for detection efficiency with the epiclccorr 
tool of SAS. We binned the light curves at 30 to 2000 s 
intervals (ensuring at least 20 events per bin) depending 
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on the source flux since it was declining. We then fit 
the fight curves to a constant, and calculated and the 
null hypothesis probability. For the Chandra CC-mode 
observations, it is difficult to correct the detection effi- 
ciency because they lack information on one of the two 
dimensions. Therefore, we assumed a uniform detection 
efficiency for the data taken with CC-mode. We then 
performed the test and the Gregory-Loredo test (us- 
ing the glvary tool of CIAO; iGregorv &: Loredolll992D 
for all the Chandra observations in Table [TJ In all cases, 
the observed light curves were consistent with being con- 
stant. 

4. DISCUSSION 

We have measured the spectral and timing properties 
of CXOU J164710.2-455216 for three years after its 2006 
outburst. We observe a clear correlation between spec- 
tral hardness and flux post-burst. In our timing analysis, 
we find a significantly smaller period derivative which im- 
plies a spin-inferred dipolar magnetic field of < 7x 10^'^ G 
(90% confidence). This is significantly lower than was 
previously reported. We find evidence that a second out- 
burst occurred between MJDs 55068 and 55832 based on 
a flux increase at MJD 55832 (Obs. 15). Next we discuss 
our flndings in relation to previous studies of this source, 
and in the context of the magnetar model. 

4.1. Correlation between Spectral Hardness and Flux 

The spectral evolution of the source after t he 2006 out- 
burst has been discussed in previous studies (jlsrael et al.l 
[2007t IWoods et al.l[20ll Based on 100-200 days of ob - 
servations. [Israel et al.r(|2007[ ) and IWoods et at] (|2011[ ) 
found no clear evidence that the spectral shape (kT and 
r) had evolved. Our analysis results for the same period 
qualitatively agree with theirs (see the first seven data 
points of the F and kT evolution plots in Figure [2]). 

However, the power-law photon index may not be a re- 
liable measure of spectral hardness because a power law 
may not be the true spectrum, especially at low energies , 
as noted by many authors (e.g. JLvutikov fc Gavrii]| l2006l : 
lOzel fc Guvct1[2U07I : [Fernandez fc Thompson 2007.) . Sev- 
eral authors showed that hardness ratio correlates with 
flux m ore strongly (e.g., lOzel &: Giiveil 120071: iZhu et al.l 
[20081) . For this reason, we plot the hardness ratio 
(^4-10 kcv/^2-4 kcv) as a function of 2-10 keV flux in 
Figure [21 The hardness ratio showed a clear correlation 
with flux even for the flrst ~100-200 days after the 2006 
outburst. This proves that the spectral shape changed 
during that period. Moreover, given the longer baseline 
of the observations, we now discern a clear trend of pho- 
ton index and kT with fiux, in agreement with spectral 
softening as the source relaxes. 

We note however that the first data point obtained 
1.7 days after the outburst epoch lies outside the trend 
for F. The source was hottest (~0.7 keV) at that 
epoch but the power-law co mponent was not hard- 
est. For a crustal event (see [Muno et al.[ [2007I). this 
may be explained as follows. When a crustal event 
occurs, the crustal fracture may implant twists in 
the magnetosphere. The twis ted fields then induce 
currents in the magnetosphere ()Thompson et al.[ [2002[ : 
[Beloborodov fc Thompso"n|[2007[ ). Perhaps the magneto- 
spheric plasma was not yet fully activated at the epoch of 



the first post-outburst observation. This would be consis- 
tent with what we infer from the change in the blackbody 
radius; the magnetospheric model predict s that the X- 
ray e mitting area decreases monotonically (jBeloborodovl 
'20101), whereas the data show an increase in the begin- 
ning (see blackbody radius plot in Figure [2|) . 

We note that the flux evolution after Obs. 15 showed 
a similar trend as during early times after the 2006 out- 
burst: the power-law spectrum hardened while the black- 
body temperature decreased and the radius increased 
from 0.67±0.03 km (Obs. 15) to 0.87±0.08 km (Obs. 16). 

4.2. Flux Evolution 

[Israel et al.l (|2007D suggested that the 2-10 keV 
absorption-corrected flux relaxation of the source fol- 
lowed a power-law decay with index of —0.28 ± 0.05 for 
approximatel y 120 days, based on Swift observations. 
[Woods et ail ()2011[ ) obtained a power-law decay with 
index of —0.306 ± 0.005 for the flux evolution for the 
first ~200 days based on Chandra, XMM-Newton and 
Suzaku observations. We find that either a double ex- 
ponential or a broken power-law function describes the 
long-term cooling trend of the source, although a double 
exponential provides a significantly better fit. For the 
first 100-200 days, the cooling trend can be described 
by a power-law decay with index of —0.22 ± 0.01. This 
is roughly c o nsiste nt wi th the previous mea surements of 
[Israel et al.[ (|2007[ ) and [Woods et al.[ (|2011[ ). Note that 
the two data points between 100 days and 200 days 
are attributed to th e second pow e r-law decay compo- 
nent in our fit while [Woods et al.l ()20lH ) used a single 
power law out to ~ 200 days. However, we find that 
the flux evolution at later times changed signiflcantly as 
is seen in Figure |TJ This behavior is obse rved in some 
magnetars' post-outburst relaxation (e.s.. IWoods et al.l 
12001 [Livingstone et al.ll20TTI : [An et al]l2012f r and can be 
qualitatively explained by crustal cooling models (e.g., 
[Lvubarskv et al.|[2002[) and /or the untwisting magnetic 
fields model (Beloborodov"2009'). 

Further, Israel et al. (2007) argued that the power-law 
component flux decayed more rapidly than the blackbody 
component did for the flrst 4 months. They argued from 
this that hot spots on the surface sustaine d heat longer 
than the external source did. However, [Woods et al.l 
(|2011l ) found no evidence that the power-law spectral 
component declined faster based on the stability of the 
power-law to blackbody flux ratio during cooling. Our 
results are consistent with the former; the power-law 
spectral component decayed faster than the blackbody 
component for the flrst ~100 days. This trend is clearly 
visible in Figure [1] Also note that the power-law de- 
cay indices we m easured are consistent with those of 
[Israel et al.[ (|2007[ ) although the flt was not good for the 
b lackbody coinponent . 

[Beloborodov[ ()2009[ ) proposed the "untwisting" mag- 
netosphe ric model to e xplain transient cooling of mag- 
netars. [Beloborodov! (l2010l) argues that the area of 
a hot spot shrinks in the "untwisting" model while it 
should not in the crustal cooling case since heat is ex- 
pected to diffuse. On the basis of observations of sev - 
eral magnetars' transient cooling, [Beloborodov[ (|2010[ ) 
argues that magnetospheric untwisting plays a signifi- 
cant role in the transient cooling. For the cooling of 
CXOU J164710.2-455216 after its 2006 outburst, there 
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is evidence that the area of the hot spot increased for the 
first < 10 days and then decreased, if we assume that the 
blackbody component in our spectral model reasonably 
represents the true thermal emission from the source. In 
this scenario, the evolution of the blackbody radius and 
the photon index imply that the relaxation is a com- 
bination of the crustal and the magnetospheric effects 
with the dominant process being the crustal cooling at 
early times. Whether or not CXOU J164710.2-455216 
is the only source that showed this behavior is not 
clear. An increase of the blackbody radius at the very 
early stages of relaxation and subsequent decrease has 
been observed in other magnetars (e.g., SGR 0501+4516 
and SGR 0418+5729; iRea et al.l [20091 lEsposito et al.l 
120101) . However, sudden hardening of the power- law 
component was not observed for SGR 0501+4516, and 
a power-law component was not detected significantly 
for SGR 0418+5729. 

With the current level of theory predictions and model 
developments, it is difficult to explain the changes of 
all the observational properties such as flux, spectrum 
and blackbody area simultaneously, and to know unam- 
biguously whether the external (magnetospheric) or the 
internal (crustal) effect dominates at any epoch. Nev- 
ertheless, here we consider a crustal cooling model in 
more detail, and show that it can reproduce the ob- 
served flux decay. In the model, we calculate the ther- 
mal relaxation of the crust after a rapid energy depo- 
sition by solving the diffus ion equatioii (see Gumming 
et al . 2012 in preparation; iScholz et al.l 120121 : lAn et ahl 
120121 for recent application s). The method of calcula - 
tion is the same as that of iBrown fc Cummind ()2009f ) 
for accreting neutron stars, but modified to include the 
effect, of strong m agnetic fields on the microphysics 
(|Aguilera. Pons fc M irallcs 2008). 

The model is essentially ID, but the effect of the 
magnetic field on the heat transport is taken into ac- 
count by assuming a dipo lar magnetic field and aver- 
aging over spherical shells (i Potekhin fc Yakovlev I [20011 : 
[Grccnstcin & Hartkc 1983). We use reasonable values 
for the neutron star mass (I.4M0) and radius (12 km). 
We then find the best set of parameters, the magnetic 
field (B) and the core temperature (Tc) and the initial 
temperature profile as a function of density to explain 
the fiux relaxation of a source (see Gumming et al. 2012 
in preparation for more details). 

We applied the model to match the light curve of 
GXOU J164710.2-455216 after the 2006 outburst. Fig- 
ure [1] shows sample cooling models that are qualitatively 
consistent with the data. The total energy and depth 
are robust, and correspond to '-^ 1 x 10'*'* ergs injected 
into most of the outer crust. The profile of the heat- 
ing must be such that the temperature profile is close to 
being independent of density, with the exact slope de- 
pending on the B that we assume. For the models with 
B = 8 X 10*3 G or B = 2 X 10** G, the temperature 
must decrease slowly with density T oc p""'^, whereas 
for i? = 2 X 10*^ G, an isothermal model reproduces the 
observed decay. 

For B = 8 X 10*3 G or B = 2 X 10** G, we find that 
the inner boundary of the heated region has to be placed 
at a density 5 x 10** g cm"^^ close to neutron drip (at 
p w 4 X 10** g cm^3) The fact that the inner boundary 
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Figure 4. Blackbody area versus 2-10 keV X-ray luminosity of 
the source after the 2006 outburst. The solid line shows Lx oc 
A^'^, where An is the blackbody area in units of 10^^ cm^, and 
dashed lines show the simple models (Lx oc A^j) in Figure 3 of 
IBeloborodo"^ 1I2OIOI '). 

is SO close to the neutron drip density may indicate that 
the heating extends into the inner crust, but that the 
temperature close to neutron drip where the neutrons 
are normal (low critical temperature) remains low due 
to the neutron heat capacity. This will be investigated 
further elsewhere (Gumming et al. 2012 in preparation). 

The fact that the temperature profile is relatively in- 
dependent of depth in the models is interesting; if the 
energy w as deposited as a fixed energy density as as- 
sumed bv lLvubarskv et all ()2002D and corresponding to, 
for example, a fixed fraction of the magnetic energy den- 
sity, the temperature profile would decrease quite steeply 
with increasing density (T cx p~°'^) and that would give a 
light curve that drops too steeply for t < 100 days, incon- 
sistent with the data. Instead, the preferred energy de- 
position profile is more like "constant energy per gram" 
(or equivalently per particle) . Grust breaking would give 
something close to this, because the maximum elastic en- 
ergy that can be stored in the crust is proportional to the 
shear modulus /i, which is proportional to the pressure 
n, so that the energy density deposited by crust breaking 
would be oc /i oc n cx p^^^ , and the temperature profile be 
T oc p^'^ , increasing only slowly with density. Our model 
profiles sit somewhere between the the fixed energy den- 
sity and the maximal crust breaking cases, suggesting 
that crustal breaking might have occurred only locally. 

Our models calculate the bolometric luminosity evolu- 
tion, but we consider only the 2-10 keV source flux since 
it is a more robust measurement than in the 0.5-10 keV, 
due to the uncertainty in A'h. If we considered the bolo- 
metric flux, the shape of the flux evolution, and thus the 
model parameters might change. It was difficult to in- 
vestigate this with the spectral models we used, since we 
are using a power-law model which has no lower-energy 
cutoff. Instead, we attempted to determine bolometric 
fluxes by utilizing the RGS model fits. 

Although the RGS model parameters were not well 
constrained (see Section 13. 2|) , we were able to measure 
the fiuxes well with this model. The RGS fluxes agreed 
very well with those of the power-law plus blackbody 
in the 2-10 keV band, and thus our crustal cooling 
models were able to flt the RGS flux evolution without 
changing any parameter in this case. The bolometric 
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Figure 5. Location of CXOU J164710.2-455216 in the P-P di- 
agram. It sits well below the B = 10^* G line. Pulsar data are 
taken from the ATNF pulsar database and magnetars are from the 
McGill online magnetar catalog. 

fluxes were higher and the discrepancy slowly increased 
ometric /F2-10 keV = 1.7 — 2.4) as the flux decreased. 
The shape of the flux evolution changed in this case, so 
one may expect the cooling model parameters to be mod- 
ified. However, the quiescent flux level also increased, 
and thus we had to adjust the constant offset of the 
crustal cooling models (core temperature). The offset 
compensated for the change of the shape, and no other 
parameters (except for the normalization which corre- 
sponds to the total energy in the model) needed to be 
modified significantly. For example, the core tempera- 
ture of the model increased by ~ 40% (1.5 x 10^ K), the 
total energy by a factor of ~ 2 (~ 4 x 10''^ ergs), and the 
inner boundary of the heated region had to be placed at 
~ 10% higher density (p = 5.5 x 10^^ g cm~^) for the 

low magnetic field model {B — 8 x 10^^ G); 

We note that the untwisting models (|Beloborodovl 
[2OO9I) can produce diverse functional forms for the fltrx 
decay and the blackbody area (Abb) evolution, and may 
as well explain the evolution of source's flux and black- 
body area during the period over which t he blackbody 
area decreased monotonically. For example. iBeloborodovl 
(|2009) shows a sample cooling curve (see his Figure 7), 
and Beloborodovj (j201 0) predicts that luminosity is 
quadratically proportional to the blackbody area, for a 
model with constant discharge voltage (see his Figure 3). 
However, we find Lx oc better describes the area 
versus luminosity relation of the source after the 2006 
outburst (see Fig. |4]). Therefore, the simple model may 
not properly describe the flux and the blackbody area 
evolution of the source. Investigating other models is 
beyond the scope of this paper. 

4.3. Timing Properties 

[Israel et al.l ()2007D and iWoods et all (|201lD measured 
the timing properties of the source with data covering 
~100-200 days after the 2006 outburst. Their measure- 
ments together with the results of this work are summa- 



rized in Table [H The measured periods at the outburst 
onset epoch agree well with one another. However, the 
period derivative we measured is significantly smaller. 
This might be due to enhanced spin-down following a pu- 
tative glitch at the outburst epoch (Israel et al.''2007), as 
has been seen in oth er sources ( IRXS J170849.0-4 00910, 
IE 2259 +586:lKaspi & Gavriil 2001 IWoods "erall l20041 . 
Indeed, IWoods etal.. (,2011. ) noted that glitch recovery 
could bias the magnitude of the spin-down rate upward 
for this source and that the spin-inferred dipolar mag- 
netic field of the source would be ~ 3.7 x 10^"^ G as- 
suming that the glitch recover y was completed by MJD 
54148. Although [Woods et al.l (|20lH ) concluded that the 
presence of a glitch could not be conclusively demon- 
strated, the small spin-down rate we measure supports 
the idea of glitch recovery. If this is the case, our re- 
sult implies that the spin-inferred dipolar magnetic field 
of the source is less than 7 x 10^'^ G, wh ich is consis- 
tent w ith B = 3.7 X 10^3 G as inferred bv IWoods et all 
(|2011l ) and comparable to those of other apparently 
low magnetic- field magnetars (e.p;., Sw ift J1822.3— 1606, 
IE 2259-H8 6; iLivingstonc et al. 2011: Rea et al.l 120121 : 
iScholz et al] |2012; Gavrnl & Kaspi 200l) and high- 
B rotation-power ed pulsars (e.g., PSR J1718— 3718: 
INg fc Kaspilf20ll . 

With this measurement, CXOU J164710.2-455216can 
be added to the growing list of neutron stars display- 
ing magnetar behavior but having relatively low mag- 
netic fields (< 10^"* G), comparable to those of high- 
B rotation-powered pulsars. Figure [5] shows the lo- 
cation of the source in the PP diagram, where pul- 
sar data are taken from the ATNF pulsar databas^HI 
(Manchester et al. 2005) and magnetars are from the 
McGill online magnetar catalog. Given similar magnetic 
fields, it appears that some show magnetar-like behavior 
while others do not. This is puzzling in the context of 
the magnetar model unless we assume higher order mul- 
tipoles or strong internal toroidal fields in some sources 
(e.g. SGR 0418+5729; IRea et al.l IMol : IGiiver et al.l 

We note that the dipolar magnetic field strengths 
measured using the standard dipolar braking relation 
may not be correct since the formula includes many as- 
sumptions such as radius, mass and magnetic inclina- 
tion angle. For example, the inferred i?-field strength 
can chan ge by a factor of two due to the inclination 
angle (.Spitkovskvl 12006 ). Also the mass (for Mns — 
1.4 ± 0.5Mq) and the radius (for i? = 12 ± 2 km) can 
change it by ^ 15% and ~ 30%, respectively. Although 
these uncertainties are fairly large, they are not likely to 
explain the larger range (three orders of magnitude) of 
inferred magnetic field strengths now observed in mag- 
netars. 

5. CONCLUSIONS 

Using archival data, we have measured the spectral 
evolution of CXOU J164710. 2-455216 in the 2-10 keV 
band for approximately 3 years since its 2006 outburst. 
We see a clear correlation between spectral hardness and 
the flux: the spectrum softens as the flux declines. Our 
timing analysis for data spanning approximately 2500 
days shows that the spin-inferred dipolar magnetic fleld 

'-^ http:/ /www. atnf.csiro.au/research/pulsar / psrcat 
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of the source is less than 7 x 10^'^ G. This is significantly 
lower than what has been previously reported and 
supports the possibility of glitch recovery following the 
2006 outburst. This result adds to the growing list 
of relatively low--B magnetars. We find evidence of a 
second outburst based on a flux increase between MJD 
55068 and 55832. Finally, fitting the flux decay with a 
crustal cooling model suggests that the cooling trend 
of CXOU Jf647f0.2-455216 after its 2006 outburst 
can be reproduced if energy was deposited in the outer 
crust and the initial temperature profile was relatively 
independent of depth. 
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